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Fermions in the harmonic potential and string theory 
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Abstract 

We explicitly derive collective field theory description for the system of fermions 
in the harmonic potential. This field theory appears to be a coupled system of free 
scalar and (modified) Liouville field. This theory should be considered as an exact 
bosonization of the system of non-relativistic fermions in the harmonic potential. Being 
surprisingly similar to the world-sheet formulation of c = 1 string theory, this theory 
has quite different physical features and it is conjectured to give space-time description 
of the string theory, dual to the fermions in the harmonic potential. A vertex operator 
in this theory is shown to be a field theoretical representation of the local fermion 
operator, thus describing a DO brane in the string language. Possible generalization of 
this result and its derivation for the case of c = 1 string theory (fermions in the inverse 
harmonic potential) is discussed. 



1 Introduction and outline 

The c = 1 string theory (see e.g. pi EJ 0] for review) is probably the simplest non-trivial 
representative in the vast hierarchy of string theories. Although it lacks many properties 
of its siblings, it is still non-trivial example, admitting interesting space-time physics. It 
offers therefore very useful and already standard laboratory for studies of general properties 
of string theories: open-closed string dualities, D-branes, tachyon condensation, as well as 
those phenomena, which are still hard to describe in higher-dimensional string theories (such 
as e.g. black hole formation). The latter problem is of the special interest, because two- 
dimensional string theory is likely to be the first model of black hole [H El E] which can be 
elaborated to such an extent, that it would address all important questions of black hole 
physics, from the initial formation to the description of final stage of Hawking evaporation 
(fully taking into account quantum backreaction) . 

Initial progress in c = 1 string theory was made in the early 90s (see reviews Ql 121 El El) 
when its "dual description" in terms of the so-called c = 1 matrix quantum mechanics 
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Figure 1: MQM in such potential has two universal regions - harmonic potential near its 
minimum and inverse harmonic near its maximum. 

(MQM) was found and the detailed correspondence between perturbative string theory and 
matrix model was developed. This allowed to use matrix model techniques for computations 
of the string theory correlators in a way alternative to the Liouville theory and provided non- 
trivial cross-check. Space-time string field theory was developed [HI IHl ^I] • The integrable 
structure of MQM allows to solve the c = 1 string theory in the large class of non-trivial 
backgrounds [HI [T21 UHl • Recently, the interest to the correspondence between c = 1 string 
theory and MQM was revitalized. Introduction of the D-branes as well as the progress in 
the boundary Liouville CFT [T^ ITB] allowed for a new understanding of the tachyon 
condensation and open-closed string duality (for a recent review see e.g. and helped 
to extend the connections between MQM and string theory to the new non-perturbative 

effects [TBI [Tnii2iaiziii22ii2nii2i!. 

Despite of all this progress and apparent simplicity of the model, the problem of black 
hole is still unsolved [^ 1^ . This may imply that some important ingredients are missing in 
usual set up of c = 1 string theory. That is one reason why it is interesting to consider some 
generalization of the above construction. In this paper we will consider such a generalization 
which can be both helpful in better understanding of the relation between MQM and c = 1 
string theory and interesting by itself. 

Singlet sector of MQM can be represented as non-relativistic fermions in the inverse 
harmonic potential. This fact can be viewed as a reason why the model is solvable. In terms 
of such fermions the most natural realization of the integrable structure of the theory can be 
formulated. Moreover, recently these fermions acquired physical interpretation [THl 12m 1^. 
Therefore they should be considered as fundamental objects of the theory rather than just 
a technical tool. So it is useful to consider generalizations of the theory which are the most 
natural from the point of view of the fermions and not only from the world-sheet string 
theory perspective. 

The simplest modification of this kind is to consider the fermions in the harmonic (rather 
than inverse harmonic) potential. This system has several interesting features. First of all, it 
has a discrete spectrum which, for example, makes much clearer the constructions of many- 
particle wave functions. This is important e.g. when one is interested in non-perturbative 
effects. Another important feature is the finite depth Fermi sea instead of infinite Dirac 
sea, which is natural for non-relativistic fermions. These features do not necessary imply 
that the theory has no dual string theory description. One can imagine the theory in, say, 
potential sketched on the Figure H Such a theory looks like inverse harmonic oscillator 
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near the maximum (apart from the fact that all states are metastable) and almost like 
harmonic oscillator near the minimum. The complete theory in such a potential can be 
considered from the point of view of the inverse harmonic potential as one of the possible 
(non-universal) regularizations. Both limits are universal and are related by the analytic 
continuation (see [211 12B] for the discussion). In the first limit the (double)scaled theory 
has infinite Fermi sea and continuous spectrum. This corresponds to the continuous limit 
of MQM from the point of view of surface triangulation and gives c = 1 string theory in 
the linear dilaton background. If one is instead interested in some other string theory, this 
feature is not necessary. Moreover, it can in principle even prevent one from describing in 
terms of the fermions a theory which in fact does not imply infinite Fermi sea for fermions 
(D-branes) whose condensate is a dual description of its background. 

Outline of the results 

In this paper we take the following approach. We try to find a bosonic theory which is 
equivalent to the system of fermions in the harmonic potential in a sense that all operators 
and their expectation values are in one to one correspondence. This correspondence is 
directly derived and can be viewed as an example of an exact bosonization of non-relativistic 
fermions.^ 

This derivation has two steps. First, one tries to bosonize fermions in terms of free scalar 
using the analogy with the theory of relativistic fermions. In the case of the inverse harmonic 
potential this gives an interacting and highly non-trivial effective field theory [HI CD] which 
simplifies only asymptotically, far from a classical turning point, where fermions are almost 
relativistic. 

Many-particle state of the fermions form a Fermi sea in the phase space of the system. The 
oscillators of this bosonic effective theory correspond to the small deformations propagating 
along the Fermi sea. These deformations are well-defined in the phase space and are two- 
dimensional in their nature. It means, in particular, that the wave function, corresponding 
to the perturbed Fermi not eigen-functions of the Hamiltonian with the harmonic 

(inverse harmonic) potential. If one tries to find a Hamiltonian for which this wave-function 
will be an eigen-state, such Hamiltonian turns out to be very complicated and non-local 

(c.f. [I3I2H1).' 

However, in both cases of harmonic and inverse harmonic potential there exist a very 
simple equivalent reformulation of the problem in terms of a system of free fermions in two 
dimensions in the strong uniform magnetic field |2H|- In such field the system is confined to 
the zero energy level of its Hamiltonian. This level is infinitely degenerate. It was shown 
in [2H1 that all wave-functions which correspond to all possible incompressible deformations 

^The exact bosonization for the case of inverse harmonic potential was first discussed in ^312111211, where 
it was reahzed in terms of quantum phase space distribution of the Fermi hquid (see for a recent results 
and applications to the study of tachyon dynamics in c = 1 string theory). This approach is somewhat 
similar to the one, undertaken in this paper. We leave detailed comparison for the future work. 

^The system is integrable and this non-local operator is related to the original Hamiltonian by the dressing 
procedure in the terminology of integrable systems. 
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of the Fermi sea belong to this leveL The equivalence with the one-dimensional problem can 
be established as follows. To consider thermodynamic properties of the system, one needs 
to choose an A^-particle ground state. According to the Fermi statistics, it is usually build 
as a Slater determinant of the first N eigen functions of the Hamiltonian. When working 
with the space of the solutions of the equation H2d = 0, to do this one needs to choose some 
ordering on this infinitely degenerate energy level. This is done by choosing any (Hermitian) 
operator commuting with the Hamiltonian and choosing ordered basis of its eigen-functions. 
If one chooses the angular momentum operator L for the ordering, the system of equations 

H2cflpE{z+, 2-) = 0, LiPe{z+, z_) = EiPe{z+, z_) (1.1) 

is exactly equivalent to the one-dimensional equation HidipE{z) = EiI)e{z) and the A^-particle 
ground state of this system coincides with undeformed Fermi sea in the phase space of one- 
dimensional system. If one chooses some other operator for the ordering, the ground state 
and, therefore, the shape of the Fermi sea will be different. But the wave-functions corre- 
sponding to all of them have one common property - they satisfy first of the equations ()1.1|) . 
That is why two dimensional representation is convenient for the study of effective field 
theory description of the system. 

As a result, instead of eigen- value problem of one-dimensional Hamiltonian of MQM one 
has a Wheeler-DeWitt type equation saying that the two-dimensional Hamiltonian is equal 
to zero. The same happens if one relates e.g. quantization of the relativistic particle in terms 
of coordinate time (after fixing some space-time splitting) to the covariant quantization of the 
same system with respect to its proper time done without any coordinate choice. This is why 
such a two dimensional reformulation can be considered to be a relativistic reformulation of 
the MQM system j2Hl- In this two-dimensional formulation different shapes of the Fermi sea 
on the MQM side all become solutions of the same two-dimensional Hamiltonian (Wheeler- 
DeWitt) equation, corresponding to the different choices of the space-time splitting. 

This reformulation allows to write down a bosonization of the system in terms of the free 
chiral two-dimensional boson. In case of inverse harmonic potential, this result is contained, 
somewhat implicitly, in jT^] and for both potentials in |2H]- 

The bosonization described above has nevertheless one unpleasant feature. The theory 
of free boson obviously has ultra-violet divergences. On the other hand, the initial theory 
of non- relativistic fermions is finite at any scale. The reason for this discrepancy is that 
the system of such fermions has the finite number of degrees of freedom. Therefore, not all 
degrees of freedom of bosonic field are independent for finite A^. To be able to reconstruct 
properly quantum-mechanical correlators for a finite, but large A^ from the field theoretical 
description, one needs to be more careful with the large A^ limit. 

The solution of this problem, described in the present paper for the case of the harmonic 
potential, gives rise to the following picture. Instead of chiral boson, we have the system of 
the full boson X plus (modified) Liouville field with the action 

5iL+x[X, 4>] = ^j d\ \dXdX + d<pd<P +J<P- /ie^'^l (1-2) 
(where h is a constant of order 1/A^). Vertex operator in this system, corresponding to the 
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local fermion operator on the quantum mechanical side, looks like 

^+ (C) = eK'^K'C)+x(C)-x{c)) 

Let us stress, that this theory is an exact bosonization of the system of fermions in the 
harmonic potential and it is derived explicitly in the present paper. In the limit when the 
Fermi step-function, Liouville field freezes to its background value and effectively 

imposes the Dirichlet boundary condition on the boson X. In this way we recover the initial 
bosonization in terms of free chiral boson. 

Surprisingly, this system looks very much like string theory. There are, however, two 
differences with the string theory in the linear dilaton background |21 C]. First, the 
exponent in the Liouville field is imaginary. This is due to compactness of the Fermi sea. 
This property can be removed by analytic continuation, which should be done simultaneously 
in the z-plane and in the space of the field 0. This continuation is tricky. For the chiral boson 
without Liouville field analytic continuation in the ^-plane relates the theories in harmonic 
and inverse harmonic potentials as well as bosonization formulas for these theories |28J. 
Generalization of this analytic continuation for the full theories is discussed in the section |S] 

Another important difference is the appearance of the background term in ()1.2|) . This 
background is not a usual curvature term as in the action of c = 1 string theory and thus 
cannot be removed from the action by a shift of field cf). Indeed, in the c = 1 string theory 
the field redefinition —>■(/) + a can be compensated by the Weyl rescaling of the world- 
sheet metric Qab Qab^^'^ ■ Thus, the curvature term can be eliminated and the Liouville 
term y/giie"^ remains unchanged. In the case involved, we do not have an analog of the ^ 

term in the action (jl.2p . As a result, the field redefinition icj) ^ icj) — would modify the 

|^|2 ^ 

Liouville potential /i e**^ — /i e**^^" (we will return to this question in Sections l2. If IHjl. Thus, 
this linear term is a real modification of the theory and it would be there also in the same 
analysis for the inverse harmonic potential. The result of this modification is important - it 
is responsible for the regular behavior of the correlators at small scales, as discussed in the 
Sections 121 13.11 This feature, which is the main reason why Liouville field is necessary in 
the bosonization of the non-relativistic fermions, is expected from the space-time picture of 
a string theory. Indeed, our analysis is the analog of the Das-Jevicki collective field theory 
approach (made for harmonic potential, in a different way and valid at any scale). This 
collective theory is believed to be indeed the space-time theory. Therefore our theory should 
have the same interpretation as well. This space-time interpretation and why our theory is 
so similar to the world-sheet theory is not yet clear for us.^ 

■^The fact that the fermions in the harmonic potential have dual string theory description was also noticed 
in [32]. When this paper was at the final stage of preparation, the work appeared, where the theory 
of "imaginary" Liouville field coupled to the free boson was claimed to be dual to the matrix model in the 
harmonic potential. In our paper we give a direct (and independent) prove of a similar fact - namely, that 
the system of free fermions in the harmonic potential is equivalent to the system H1.2|l which is almost (but 
not completely!) the same as the one of |3S1- One important difference between our theory and that of 
is that in the latter background term is treated as a usual curved metric background in Liouville theory. 
Thus it can be removed by a field redefinition. In our case it is not true, as it was discussed above as well 
as in the end of Section ITTl The importance of this modification is that it makes correlators finite at small 
scales. This makes two theories physically different. Due to the same reason, "duality map" seems to be 
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Figure 2: Typical density distribution as a smoothed out step-function. The plot is shown 
for = 64 particles. 



2 Collective field theory for the fermion density 

We study one-dimensional fermions in the harmonic potential. As it was discussed in j2H] 
and reviewed in Section ^ this theory can be equivalently reformulated as the system of two- 
dimensional fermions in the magnetic field. In this section we will show that it is equivalent 
to the "imaginary" Liouville field theory. This field theory is in its spirit a collective field 
theory [311 |H], written in terms of the density of fermions, rather than describing the dynamics 
of individual fermions. 

We begin with the following family of ground states of the system of fermions:^ 

1 ^ -^El^f 
^o(^i, ■■■.zm) = {zu ...,zn\0) = —=- l[{z, - z^fe " ' (2.1) 

The discussion below will be concerned with the statistically large number of particles. The 
thermodynamic limit of this system is defined as — > oo. It is well known |^|37j that in 
this limit ground state (|2.1|) describes fermions forming round droplet of radius with the 
constant density equal ^ inside and zero outside. Namely, the average of density operator 

different as weff. For example, local fermion operators correspond to the vertex operators on the "string 
theory" side in our case. The two-point function of vertex operators is non-singular, in agreement with 
the quantum mechanical result. At the same time if one calculates two point function of the same vertex 
operators in the model suggested by '33', it will have an ultra-violet singularity usual for a conformal field 
theory. It is therefore unclear how to build the operators representing fermions correctly in the model .33, . 
Unfortunately, in the space-time interpretation of the theory is not discussed. If the proposed there 
duality is correct, it gives a world-sheet rather than a space-time theory. 

"^We introduces an additional parameter (3 = 1,3,5,... In case (3=1 (free fermions in the magnetic 
field) the wave function p. 1(1 is a Slater determinant of the first N monomials 4'o ~ det(z*^^), which is 
just the ground state of the system of N non-interacting fermions, ordered with respect to their angular 
momentum. Case /3 > 1 corresponds to a more complicated interacting system: it is a celebrated Laughlin's 
wave-function jj^, describing fractional quantum Hall state with the filling fraction /3~^. Naively, case /3 > 1 
does not correspond to the singlet sector of MQM and thus has no immediate connection to our problem. 
However, surprisingly, all the derivations work exactly the same for any /3, so we can keep it "at no cost" . 
Moreover, as we will see later (Section O, (3 plays the role of the radius of compactification for one of the 
fields in the collective theory of the system (|2.1|) . 

It is interesting to note that the theory with /3 > 1 can arise in the non-singlet sectors of the MQM |36j . 
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p{\z\) = {0\p{z,z) |0) has the profile, shown on Fig. |21on the preceding page. Notice that 
the width of the "edge" (transition region between p = const and p = 0) is proportional 
to 1 in our units (and is small compared to the radius). In many cases it is enough to 
approximate this profile by the step-function: p{\z\) = ^9{\/N — \z\), however, one should 
take into account, that this is a non-uniform approximation: for \z — ^/N\ -C 1 difference 
between the real distribution and the step-function is always of order one. 

Wave-function ()2.1|1 is not normalized. The norm of the ground state 1)2.11) is defined by 

tn = J (fzi. . .(fzNl'^ol'^ (2.2) 

ttv can be interpreted as a partition sum of the Coulomb plasma, by noticing that eq. ()2.2j) 
can be written as 

Zn{P) = j (fzi . . . rf^z^l^ol' = j d^zi . . . d^N e-^^^(^i'--^^) (2.3) 

where 

N 

Sc{zu...,zn) = - J]21n|z, -2:^1 + (2.4) 

i<j i=l 

is the classical energy of a 2D Coulomb plasma of particles of charge = 2 moving in the 
uniform background of particle density po = -. This Coulomb plasma analogy allows for a 
simple and transparent way to understand what kind of collective field theory one should 
expect. First of all, it is easy to show that the density distribution has be approximated by 
the step-function (see e.g. [33 EZl EH EH] in the context of quantum Hall effect or [101 HI] in 
the context of matrix models). For many problems step- function-like density is enough, but 
let us stress once again that it is a non- uniform approximation on the edge of the droplet. 

This approximation can be significantly improved by noticing that ()2.3p looks like a 
plasma with the temperature Indeed, let us introduce the electric potential ip{z,z), 

satisfying the 2D Poisson equation 

A(p = 47rep{z,z) (2.5) 

where p = Pp — pQ is the total charge density created by both particles and background in the 
system. The key approximation, often used for such a system is to say, that the distribution 
of the mobile charge in the external electric field is given by the Boltzmann law (recall that 
the temperature is 

Pp{z,z)=pe-^'^^^^'^ (2.6) 

where p is the fugacity. Combining equations ()2.6|) and ()2.5|) we find that the scalar potential 
(fi satisfies the equation 

Aip = 47re (/ie-^"^("'"") - po) (2.7) 

For our purposes the parameter p is to be chosen so that the total number of particles be 
equal to 

p I d^z e-^'='^("'"") = N (2.8) 
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Equations ()2.7p and ()2.8|) constitute the closed set of equations which determine the density 
distribution of the fermionic droplet in our approximation.^ Solution to this equation gives 
rise to a smooth function, which serves as a better approximation to the real density, com- 
puted over ()2.1|) . Loop expansion of the field theory (discussed in the Section H)) will give a 
uniform approximation. This smoothing out is responsible for getting rid of the ultra-violet 
singularities in the fermion correlators (c.f. Section IXTj) . 

Eq. ()2.7p is clearly of the Liouville type (modulo the constant term). Thus, we can expect, 
that written in terms of field (p (which is, essentially, the logarithm of the fermion density), 
one should get the collective field theory, related to that of Liouville. 



2.1 The Liouville field 

To derive a field theoretical reformulation of the quantum mechanics of fermions ()2.1|) . 
we will show that in the large N limit the grand partition function of the system ()2.4j) (or, 
equivalently, ()2.1|) ) can be calculated by means of field-theoretical path integral.® The saddle 
point equation of this theory will reproduce (j2.7|) . 

We define grand canonical ensemble of the Coulomb plasma via^ 

Ar=0 ■ N=0 ' N=0 ' 

The first step towards the field-theoretical reformulation is to rewrite the Hamiltonian ()2.4p 
in terms of a correlator in a free boson field theory: 

/ N \ 

\^=' I So 

where the average 

(•••)5o = ^ /©^...e-^oW, Z^={l)so (2.11) 



is taken with respect to the Euclidean Gaussian action of a free massless boson' 



5oM = ^ / d'zd(l>d(l> (2.12) 



^This equation can be derived from the Ward identities in the context of normal and complex matrix 
models 01], which is, of course, closely related to this problem. 

^Such approach was first used in 42. to obtain a diagram technique for a description of phase transitions 
in various systems. 

Us tn ~ for iV ^ oo, formal sum (|2.9() is divergent. This, however, can be cured by an IR 

regularization. For details, see Appendix 1X1 and work f43j . 

^One has to be careful with the presence of both UV and IR divergences in this action. We will ignore 
all these subtleties for now, referring reader to the appendix and to the future work (43) . 



8 



with both short and long-distance cutoffs imphed (see Appendix^. For the normahzation 
T/v (eq. (j2.2p ) this representation gives 



TN= f dh,...d'zN/l[e''f'^'^''^^-''\'A (2.13) 

\i=l / So 



Consider now the grand partition function ()2.9|1 . Substituting the representation ()2.13j) in 
it, we find 



1 . , r [ r 



N=0 

The sum over N is easily performed yielding 



TV 



(2.14) 



where 



(2.16) 



Equations ()2.15|) and ()2.16|) define the path integral representation for Z. The coupling 
constant 

/i = 2</5 (2.17) 
has the dimension of inverse length squared.^ 

We have shown that the grand canonical ensemble of fermions can be rewritten as the 
field theory ()2.16|) . Several comments are of order here. First, we have derived the theory, 
similar to the usual Liouville field theory, however with an imaginary (rather than real) 
exponential term - Liouville theory with imaginary coupling. For the sake of shortness we 
will often call it imaginary Liouville theory and sometimes simply Liouville theory, whenever 
it does not cause a confusion. 

As already discussed in Introduction, this action has an apparent distinction from the 
ordinary Liouville theory. Namely there is a background term e~^'^' . Notice, that by the 
field redefinition — > — i(3\z\'^ this term can be traded for the linear term in the action 
of the form 2i(3(f) (as in (|1.2|) ). which resembles a background curvature term in the usual 
Liouville action: 

Sim = ^j(fz^ {dct>dct) + R{g)ct) + /ie-^) (2.18) 

Yet, it is not a curvature (unlike that in ()2.18|) ). Indeed, in the action ()2.18|1 one can 
perform the shift of variable — > + 0", such that R — 2V^o" = 0. Then the curvature 
term would disappear, while Liouville terms would remain the same, as an additional factor 
e'^ would be absorbed into the Weyl rescaling of the metric gab — ^ Qab^^'^ ■ On the other 
hand, in the imaginary Liouville action ()2.16|) there is no analog of the Weyl symmetry 
and thus the background term constitutes the physical modification of the action ()2.16|1 . 
This term is responsible for making short-distance behavior of the correlators, computed 

^/x is a running coupling constant satisfying the RG equation 4^^f^ = P. See Appendix El for details. 



9 



in the theory ()2.16p finite. We will discuss the importance of this term, as well as other 
(dis)similarities with the usual c = 1 string theory in Sectional 

Notice also, that the problem, mentioned in the footnote [7| on page|Hl translates in the 
language of the theory into the well-known fact that one does not have an analytic 

expansion in parameter (as could be erroneously concluded from the definition ()2.9|) ). 

At last, it can be shown (if one takes care of all the subtleties, which arose in the derivation 
of Sect ion im c.f. AppendixIXj) the correlators computed with respect to the grand canonical 
and canonical ensembles with the number of particles given by iV = {N)gce will be the same. 



2.2 Operators in the imaginary Liouville theory 

Field theory (j2.16|) should be an equivalent description of the original system of fermions 
in the magnetic field with the ground state Thus, we need to identify the vertex 

operators in this theory, whose correlators reproduce the quantum-mechanical averages of 
the fermions. 

We begin with the fermion density operator, given in the first-quantized formalism by 
{^\p{C)\^) = j d'z,...d^ZN\'^{zi,...,ZN,0\' (2.19) 
On the ground state this operator acts as 

N 

(0|p(C)|0) = e-^l^l' / d'z,...d'z^ll\C-z,n^o{z^,...,z^)\' (2.20) 

i=i 

To find a path integral representation of the density operator ()2.19|) we need to pass to the 
grand canonical ensemble in the correlator (0| : p{zi) . . . p{zn)' |0): 

1 °° ]v\ 1 r 

(01 : p(.0 . . . p(.,): |0> = ^ E W.W^-^ J ■ ■ ■ 

Substituting the exponential from 1)2.211) with eq. ()2.10|1 and performing the sum we obtain 
the n-point density matrix in terms of the average in imaginary Liouville theory: 

(0|:p(^i)...p(^„): \0) = C /f[e''^^^-'^^-^\^A (2.22) 

where the average is taken in the theory with the action ()2.16|) . Comparing the right hand 
side and the left hand side of eq. ()2.22|) we identify the density operator in the Liouville 
theory: 

p(z, z) = Ji^e*<^("'"")-^l^l' (2.23) 
2tcp 

where p is defined in eq. ()2.17p . 
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Thus we can reproduce in the field theoretical formulation density-density correlation 
functions of the quantum mechanics of fermions. In particular, saddle point equation of (j2.16|) 
(which coincides with (12. 7|) as expected) reproduces qualitatively the correct behavior of 
the density. Indeed, it can be shown (see |43j) that it looks as a step function smoothed 
out at the edge (see qualitative picture on the fig. El on pageEj). Other properties of the 
density correlators will be discussed later. Here we would like to stress that to have the full 
correspondence between the quantum mechanics and Liouville theory we are still lacking at 
least one important ingredient. Namely, an important object in quantum mechanics of the 
system of many fermions is a local fermion operator. For example, instead of the density 
p{z,z) we could calculate fermionic two point function {ip^{z)ilj{z'))}^ Important feature of 
the fermion operator is that it is chiral. In the first quantized formalism the operator V'^(C) 
multiplies the wave function by {zi — C)-^^ 

N 

^qm(C)^o(2;i ...zn) = ^0(2:1 ...znX) = - 0^0(2:1 ...zn) (2.24) 

i=l 

The only ingredient at our disposal so far - the Liouville field 0(2;, z) - is not chiral: insertion 
of density operator exp(i0) corresponds to the multiplication of the ground state by the non- 
chiral object Y\ \zi — CP (c.f. ()2.20|) ) and the solution of the saddle point equation if (^o(k|). 

i 

Thus, it turns out that to bosonize the fermion operators one needs to have an additional 
chiral field X{z) as it was discussed e.g. in [121 EHl- In the following sections 131- 14. II we will 
remind this bosonization procedure starting from the quantum mechanics of the fermions 
and will show how to blend it with the construction of the section ITTl In short, the answer 
will be the following. We will introduce an additional free scalar field X{z, z) in the Liouville 
field theory: 

5iL[X,0] = 1-Jdh \^dXdX + d(Pd(P - /ie'^^l^l'} (2.25) 
In terms of the fields (f){z, z) and X(z, z) the local fermion operator in the theory ()2.25p will 



,4(0 . exp I M«l±i^Kl^ 1 ,2.26) 



be defined as 



where X{z) and X{z) are chiral and anti-chiral components of the field X{z, z). 

Expression (j2.26j) can be understood as follows. As we have seen in (I2.2(J|) e*"^ multiplies 
the wave function by the factors H IC ~ ^jP- Then, using this fact and eq. ()2.24|) (together 
with its anti-holomorphic counter-part), we can see that the vertex operator ()2.26|) turns the 

^°0f course {if)^ {z)'ip{z')) — > p{z, z) as z' z. 

^^We will omit /3 in the rest of this section to make presentation more transparent. Full derivation will be 
presented in the section UTI 
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wave function \l/o(-2i; • • • ; Z]^) into 



V^+(C)^o(^i,...,^iv) = (^ni^^-^l'^j ^o(;^i,...,;^7v) (2.27) 

N 

= J](C-^.)^o(^i,...,^iv) (2.28) 

i=l 

that is is indeed local fermion operator! 

To summarize, we see that the theory (|2.25p is equivalent to the quantum mechanics of 
N non-relativistic fermions in the uniform magnetic field (including the correct structure 
of the vertex operators). The relation between two theories is very similar to the usual 
bosonization. Indeed, we know that for relativistic two-dimensional fermions there is an 
exact bosonization. It was shown in IT3] that the infinite number of the non-relativistic 
fermions in the harmonic (and inverse harmonic) potentials can be described in terms of free 
chiral two-dimensional boson. Nevertheless, for any large but finite number of particles N 
such bosonization is valid in the infra-red regime only. The reason for that is quite simple. 
The theory of 2d boson is a relativistic field theory and it has ultra-violet singularities at 
small scales. These singularities are of course absent in the quantum mechanical system of 
the finite number of non-relativistic fermions. Thus, to have the exact bosonization valid 
at all scales, we need more complicated theory which is finite in the ultra-violet. In case at 
hand this appears to be the theory ()2.25p - ()2.26|) . Notice, that in spite of the similarities 
between the actions ()2.25p and that of c = 1 string theory, they are quite different in nature. 
The main difference is the background term e~^'^' in the action (j2.25|) . which (as discussed 
in Section f2.2l after eq. (j2.28|) ) and is responsible for UV regularization of the theory. As 
discussed previously, this term cannot be removed by the field redefinition and thus represents 
the modification of the theory and not just a curved background of Liouville theory. 



3 Naive infinite bosonization 

In order to derive the action (j2.25j) - (j2.26p we first review here the derivation of the bosoniza- 
tion procedure for the infinite number of non-relativistic fermions in the harmonic potential 
in terms of the free chiral boson ISH] • 

We begin with the system of fermions, described by the following ground state: 

1 ^ 

^o{zi,...,zn) = {zi,..., zn\^) = Ylizj - Zi)" (3.1) 

^^This bosonization is different from the standard Das-Jevicki theory because it is not asymptotic in 
the phase space, but is vahd everywhere and is reahzed in terms of two dimensional free boson. 

^•^Because aU the wave-functions of the LLL are holomorphic (apart from the common factor e"^'^'^'^), we 
will be omitting this factor in the rest of the paper, working in what is known as holomorphic representation. 
Thus eq. is a holomorphic representation of the ground state H2.1|l . 
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Figure 3: Simplest examples of analytic curves with non-trivial t^s. On the left: ellipse 
(only ^2 7^ 0) and on the right - curve with the Z3 symmetry {t^ 7^ 0). 



As discussed in Section ITT] after eq. ()2.4p . the particle density distribution, described by this 
function, has the typical shape presented on the figure |21 on page IHl In this section we are 
going to approximate this density by the step-function so that the state ()3.H) describes the 
circular droplet of the radius \/N. The low-lying excitations of the system ()3.H1 are incom- 
pressible deformations (i.e. those, preserving the total area of the droplet, while changing 
its shape). 

To construct a field theoretical description of these degrees of freedom, we will implement 
the following strategy [2HI EH]: first we introduce the classical sources, describing area- 
preserving deformations of the system ()3.1|) . Then we will quantize them and construct their 
Hilbert space TC. 

Classically, the following set of wave-functions describes any incompressible deformation 
of the state (jSH): 

N 

{zu . . . , zM = ^0(^1, ,zm)1[ e"-^^^^ (3.2) 

1=1 

where u^z) = YlT=i^kz'' is an entire function. Thus any state \tk) belongs to the LLL. It 
is easy to understand why these functions describe incompressible deformations. Density 
average, computed over the state is constant inside the contour 7, that has the same 
area as the ground state (|3.1|) . and whose shape is parameterized by the set of the so called 
harmonic moments (c.f. [38, 39j): 

tk = -\ I h = {hf; k>o (3.3) 

irk J z" 

ext 7 

These moments have simple geometric interpretation. Circle with the center at the origin 
has all harmonic moments equal to zero. Examples of the simple shapes with non-zero 
are given on the figure El on this page. Under certain assumptions t^s serve as coordinates 
in the space of analytic curves jUmOj. Thus, at least quasi-classically states ()3.2|) clearly 
describe all possible deformations of the droplet. 
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On the quantum level, to construct the Hilbert space Ti of incompressible deformations^'* 
of the system (j3.ip it is convenient to define the so-called t -representation (we refer to PHI IHU] 
for details) - isomorphism of Hilbert space Ti. and that of the free chiral boson. 

To construct this isomorphism, we consider the matrix element between an arbitrary 
state 1$) G Ti and \tk), which we will call i -representation of the state |$) 

N 

$(4) = (41$ ) = / d?z^ . . . d^ZM n e^^^'^^^oizu . . . , %)$(^i, . . . , ^7v) (3.4) 
The scalar product in t-representation is defined as 



(3.5) 



t,.=0 



It can be easily shown |2H1 1^^] that for any two states |$i) , |$2) G 'H their scalar product 
in ^-representation (i.e. / d'^zi . . . d'^Z]\f^i{zi, . . . , Z]\f)^2{zi, . . . , zjy)) is equal to that defined 
via (|3.5|) . Thus t-representation defines an isomorphism of the Hilbert spaces. The t-repre- 
sentation of the ground state (|3.1|) is 

N 

^o(4)= / d^zi...d^ZN ne^"^''^l*o(^i,...,^iv)|' = rjv(0,4) (3.6) 

i=i 

In case P = I function Tjq{xk,yk) appearing in ()3.6p is a tau-function of the 2D Toda lattice 
hierarchy , while for /3 > 1 it lacks interpretation in terms of integrable systems. However, 
in the large N limit this function is given by |^ IHS} 1^ 

r^(4,4) = e^^(*'='*'=) (3.7) 

where -F(tfc, 4) does not depend on /? and is a logarithm of the tau-function of dispersionless 
Toda lattice hierarchy (c.f. j^EHlIll])- Thus, we see that the integrability of the system 
is recovered in the large N limit even for (5 > 1\ An important property of the function 
F{tk,tk) is 

d'^F(tk,tk) 



tk=0 
tk=0 



dtn dtr, 

with all other first and second derivatives vanishing. 

It is easy to define the structure of Fock space in Ti in t-representation. Define operators 

Id 

a^k = ^/Pktk, Ofc = ^=-^, [ok, a-n] = kSk,n, k>0 (3.9) 

VP (ytk 

From the properties (j3.8|) of the function F{tk,tk) it follows that 

a^\Q) ^ ^lEplhl. =0 vn>0 (3.10) 

tk=0 



dtn 



^^A\l the deformations we consider in this paper are incompressible, therefore we wiU often caU them 
simply "deformations of the ground state" . 
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i.e. state ()3.6|) is a Fock vacuum with respect to the operators ()3.9|) . States a_fc |0) correspond 
to the quantized version of elementary deformations with the harmonic moment (as those 
on the figure El on pagell3|). In the space 7i there exists a basis of symmetric polynomials 



(3.11) 



tk=0 



N 

For example, state |0, . . . , 0, rifc, . . .) corresponds to (^ 2^)"'=\E'o(-2i, • • • , zn)- Property ()3.8|1 



means that the states |n) can be represented as 



In) = . |0) (3.12) 



A function uj{z), entering the definition ()3.2j) of an arbitrary incompressibly deformed state 
\tk), can be expanded in the Taylor series 

n 

Combining ()3.9j) . (j3.1()j) . ()3.12j) and (j3.13|l we see that Hilbert space H has the structure of 
Fock space of the free chiral hoson-}^ 



k>0 



iflologC + a^o (3-14) 



Let us discuss zero modes Oq and Xq of the operator X{(). We know that in our theory there 
should exist a set of vertex operators Vp{(), which in the first-quantized formalism act as 
follows: 

N 

VpiOM^i, ...,zn)=1[{C- zjf<ilo{zi, ...,zn), P=l,...,(3 (3.15) 

i=i 

In particular we are interested in the operator Vp = ip'^iC) - local fermion operator, that of 
adding a particle at a point (. It is natural to expect this operator to be realized as 

V;+(C) = :e*v^^(^): (3.16) 

This will indeed be so, if we introduce a zero mode: qq = y/j3N and its conjugated Xq: 
[oq, Xq] = 1. Note, that the fact that is quantized means that the chiral boson is compact- 
ified: X = X + 27cR. Its radius of compactification can be determined from the quantization 
condition for ao to he R = ^ (indeed ^(Ce^^^) = X{C) + 27ry^N = X(C) + 27rRN). 

^^Notice, that although for /3 > 1 ground state corresponds to some interacting system of fermions, 
the theory of its low- lying excitations is still the /ree! 



15 



Let us return to the local fermion operator. Notice, that we can write 

^+1 (0| ^HO |0)^ = ('""rr, (^0, - (3.17) 

where ket and bra states in the left hand side should be understood in N and + 1 par- 
ticle sector correspondingly. Eq. ()3.16|) gives an exact bosonization for infinite number of 
fermions. By analytic continuation one can thus show that the c = 1 MQM also admits 
exact bosonization in this limit (c.f. [T3| I28j). 

Finally, from the definition of the scalar product ()3.5|) it follows that 

{ak^ = N'' a_k] = ao, = Xq (3.18) 
and thus we can write from the conjugated version of ()3.14j) : 

^(0= k"'(C)l^ = :e-^^''(^): (3.19) 



for the local fermion operator of the annihilation of the fermion at the point Notice, 
that the expression ^ is 

the radius Vn (c.f. 1471 



that the expression ^ is the point, symmetric to the point ( with respect to the circle with 



3.1 Density matrix and UV divergences 

It is instructive to compute a density matrix to demonstrate the limits of applicability of 
the approach of the previous section. The density matrix (or fermion-fermion correlator) is 
defined as 

p{z, w) = (0| ^+(z)V'(w) |0) = (0| : e'v^^(^): : e-'^""^^^: |0) (3.20) 

In computing ()3.2()|1 in t-representation, we would have to calculate quite non-trivial scalar 
product between the states (tkltpiw) |0) and {ti^lip^z) |0): 

(01 10) = (3.21) 

(factor {zw)^'^ comes from the zero mode contribution). Functions r^v are functions of 
infinitely many variables, which can be defined only implicitly. Therefore, computation of 
expression (|3.2H) may seem to be difficult. However, it is possible to use the following fact. 
Recall |1H] that one can interpret the right hand side, of eq. ()3.7p as the Dirichlet boundary 
state on the circle of radius VN : 



r7v(4,4)=exp(/3j]A;iVVfc + ---) = (^.22) 

where the last equality is understood in terms of operators ()3.9j) . ()3.18j) and their anti- 
holomorphic counterparts (another set of bosonic modes , realized as in ()3.9j) with respect 

^^Notice, that strictly speaking the second equaUty in eq. H3.19|l is not true - it contains some additional 
c-number, coming from the correct treatment of the zero modes. In this case there this factor is ^ 
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Figure 4: Profile of tlie fermion correlator in the radial direction. Function sharply peaks at 
r = 1 and falls off very fast (the case of (3N = 100 is shown). 



to another variable tk'. a^k = ktk, etc.) acting on the "vacuum" (constant function). Then 
we can rewrite (j3.21|) as 

p{z,w) = :e^v^(^(^")+^(-)): = ; e*^/^(^(^")+^H): e^^(*-*^) (3.23) 

From this representation it is trivial to compute that^^ 

p(.=,.) = e-l«H'^N'._J^2l!^ (3.24) 

To understand the meaning of this correlator, we can now put z = y/Nre''^'^ and w = 
\/iVre*^^. Then correlator ()3.24|) becomes 

-(3N{r^-logr^) 

/^(^i-^^)= (l_e-^(^.-^.)). ^3.25) 

r-dependent part of (j3.25j) peaks sharply around r = 1 and goes to zero both as r — * and 
r ^ oo (c.f. Figure E]on the current page). On the other hand, while \6i — ^2! ^ 27r, one 
can rewrite the ^-part as simply jg^zj^- Thus ()3.25p looks like a delta-function on the edge 
of the droplet times the two-point correlator of field theory of the chiral scalar. 

The immediate problem with this correlator is that it is divergent on the surface of the 
Fermi-sea (i.e. when \z\ = \w\ = a/ZV or r = 1) as 2; — w. We know that this is not the 
case in the original theory of fermions, as the answer ()2.20j) for {p{z, z)) and it is finite. The 
reason for this bogus divergence has already been discussed throughout the paper: we tried 
to bosonize the system of non-relativistic fermions with the relativistic scalar fields. For this 
we need strictly infinite number of fermions. For any arbitrary large but finite number of 
fermions the procedure, discussed in the sections 0V l3.ll breaks down, as basis states ()3.1H) 
become linearly dependent (indeed, for any finite there are only finite number of symmetric 
polynomials of variables). The other way to see where the approximation breaks down is 
to see that the quadratic approximation of tau-function ()3.22|) fails for index k ~ (we 
refer to Appendix^. Thus, when points z and w move close enough (distance of the order 
l/\fN) to each other, the simple theory of the scalar field does not suffice and one needs 

^^In eq. H3.24|l we restored factors e"'^'^'^ and took into account additional factor, mentioned in the 
footnote El on the preceding page. 
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to embed it into more general (UV- finite) theory. As we have aheady argued in the end 
of Section |21 correct short- distance limit for the density correlators can be achieved on the 
field-theoretical side by using (modified) Liouville theory. Below we are going to combine 
two approaches (that of the Section |21 with that of the current section) to construct the 
correct fermion operators and derive the theory ()2.25|) - ()2.26p . 



4 Imaginary Liouville theory plus free boson 



Our strategy during the next steps of the derivation of the eqs. (j2.25j) - (j2.26|) will be to 
apply the construction of Section |21 using for the densities and the normalization r^v the 
field theoretical representation in terms of Liouville field found in Section |21 To do so we 
need to generalize this representation for the Ti\f{tk,ik), i.e. for the case of an arbitrary 
incompressibly deformed background \tk): 



{zi, . . . , \tk) 



N 

II' 



/3 E 



(4.1) 



with io{z) as in ()3.2|) . The norm of this state is written as 



(fzi . . . (fz^e 



where 



-2 5^ log 



Zj z. 



i<j 



N 



\z. 



(4.2) 



(4.3) 



Going to the grand canonical ensemble, we arrive to the following modification of the theory 



^ 1 f 

Z{tk,ik,fJ') = ^ j^^{tk\tk)N = Y' J '^(^(^ 



where 



'S'iL [0, tk, tk] 



N=0 
1 



27^(3 



d'z 



H e 



i0+/3[a;(2)+il.(5)]-/3|^|2 



(4.4) 



(4.5) 



with /i is given by the same expression ()2.17p . It is convenient to use a slightly different 
representation of the same action. Consider the following shift 



{z, z) — > 0(z, z) + il3[uj{z) + uj{z)] 



Then the action becomes 

•S'iL [4>, tk, tk] 



yU e 



i<t,-/3\z\^^ 



(4.6) 



(4.7) 



^^This idea has been first emphasized in works [231201021) where the UV-finiteness has been achieved 
essentially by using non-commutative field theory. 
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However, the shift ()4.6|) changes the boundary conditions for the field 0. New field does not 
vanish at z = oo, but rather has there a singularity (that of the entire functions uj{z) +Cd{z)). 
These boundary conditions can be incorporated in a form of the boundary term 



St 



An 



dz 6 duj — (b dzd) dcu 



dz cudcu — (b dz ujduj 



(4.8) 



where the integrals are taken around z = oo. 



4.1 Coupling with the free boson 



As a starting point, notice that the boundary terms in ()4.8p look like classical sources (given 
by duj{z) and du}{z)) coupled to the field (j)[z,z). In the logic of the sections OHSUl we may 
try to "quantize" these sources, by promoting them to the creation operators of some field. 
Essentially, we have already done that in eq. ()3.9|) . but we have two sets of variables {tk and 
tk) in eqs. ()4.5p - ()4.8|) . which means that the whole scalar field X{z^ z) should appear. 

To begin with, we compute average of the vertex operator of the Liouville field: 



z 



a - a 



(4.9) 



where a is an arbitrary number and Z is defined by ()4.4|) . Careful treatment of the grand 
canonical ensemble (GCE) in our case shows that all the correlators, computed in the GCE 
and in the micro-canonical ensemble with the = {N)gce coincide (as N — > oo) (for details 
see Appendix |X1 and jlS]). Thus, we can simply substitute the right hand side of eq. (j4.9|) 
with Tjyitk — -^,ik — ■^), understanding A^ as a function of chemical potential /i in a sense 
just described. Namely, if we introduce two sets of Heisenberg operators (as in Section EUl):^^ 



a-k 



a-k 



^/Pk tk, 



ak 



(3k tk, ak 



J_d_ 

T^dfk 

1 d 



k6, 



k,n , 



k>0 



[ttk, a-r, 

[ak, a_n] = k6k,n, [0'k,0,n] 



(4.10) 



VPdtk' 

define the corresponding scalar fields (compare with eq. ()3.14|) ): 



^(0 



fc>0 
k>0 



ak 

ke 

dk 

ke 



C^a-k 
k 



+ zero modes 



+ zero modes 



(4.11) 



and denote by |0, 0) the Fock vacuum of ()4.10|) : |0, 0) = 0, |0, 0) = for A; > 0, then one 
naturally defines a t-t-representation (as a straightforward generalization of construction of 

""^^ Throughout this section we have chosen to ignore zero mode contributions. This is done to keep the 
presentation as simple as possible and avoid proliferation of many non-essential factors. It should be straight- 
forward to restore all the zero modes contributions. 
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Sectionini). Now we may interpret r/v(tfc — -^,tk — ^) as a state \tn) in the Hilbert space 
Ti. in t-t-representation: 



, - I , a - a , 



— 01 Ol 

(tfc,4|r^(a^-^,a^-^)|0,0) (4.12) 



Define the the vertex operator 

In t-t-representation it acts on the function (j4.12j) by shift in both and 1^: 

(3 a 13 



{Q,Q\^|J+{CX)\rN) = TN 



a 



2kC^ kC^' 2kC^ kC 



(4.13) 



(4.14) 



If one chooses a = f , then eq. ()4.14|) is equivalent to 



;v+i (0|V^+(C) |0) 



N 



(4.15) 



(fzi . . . (fzN'^o{zi, . . . , ZN)tp^ iO'i'oiZi, ...,Zn) 



i.e. eq. fj4.15p indeed describes the average of the local fermion operator. 

Combining the results of ()4.9|) . ()4.12|) . ()4.13|) and ()4.15|) . we arrive to the conclusion that 
the vertex operator 



(4.16) 



is the local fermion operator in imaginary Liouville theory, coupled to the free scalar: 



^iL+x[X, <P] = ^j d'^z \dXdX + i 



^iJ4>-f^\z? 



(4.17) 



where the only interaction between the fields (f){z^z) and X{z^z) is through the boundary 
terms, as in ()4.8|) : 



Ait 



<^ dzcp dX 






8m 



dzXdX 



dzXdX 



(4.18) 



Several comments are in order here. First, we could have written the same expression 
entirely in the path integral formalism, where it would look like 



^+((,C)) = / Me-^«[^1 . . . e' 



;v^[x(c)-x(c)] 



V(j) e 



-5iL[0A,ffc]g§<A{C,C) 



(4.19) 



(where S'o[X] is an action of free massless scalar field X and tk, tk are the positive Fourier 
modes of the field X). In the functional integral over (p positive frequency modes of the field 
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X enter only via the boundary terms ()4.8|) . One can perform the functional integration over 
to obtain a tau-function Tj:q{tk-,ik). The remaining integral over X with such insertion is 
equivalent to the operator computations (j4.14j) - (j4.15|) . 

One can understand the meaning of this insertion as follows. In operator formalism tau- 
function ttv is an element of the Fock space of the field X, created by Tjq^al^al) |0,0). If 
one uses the naive X — > oo limit and keeps only quadratic in part of r^v, one can regard 
state Ittv) ()4.12|) as a Dirichlet boundary state on the circle with the radius a//V |i5] : 

r7v(4,4) =exp(/3^fciV'=tfc4 + ...) =exp(j]^) |0,0) (4.20) 
(compare with Section ing) . Then eq. ()4.15|1 becomes 

(0, 0| . . . et^[^(f)-^K)] |D)) = (0| . . . e^v^^^'^) |0) (4.21) 

In this approximation we recover bosonization of Section El and the theory ()4.17|) reduces to 
that of the free chiral scalar in the boundary state formalism. 

5 Comparison with the c = 1 string theory 

As it was discussed in the introduction, although the theory we derived resembles very much 
the usual c = 1 string theory, it has very different properties and interpretation. First of 
all, instead of Liouville field it contains that of with imaginary charge. One could argue 
that if we repeat this derivation for the inverse harmonic potential, we could get usual 
Liouville field. To do this we need to continue analytically from to icf). Naively if one does 
such a continuation, the kinetic term in the action changes its sign and the action becomes 
unbounded from below. However, as it was argued in to relate harmonic and inverse 
harmonic potential cases, one should make also analytic continuation in z plane to get real 
variables 2+ and Z- instead of complex (conjugated) variables z and z. Then the action 
becomes Lorentzian and the analytic continuation in </> does change its properties. Another 
problem is related to the behavior of the potential. If we consider its dependence on the 
constant 0, in the complex plane there is a line where the potential is unbounded [SD] • This 
may suggest that two theories are not equivalent and the derivation for the inverse harmonic 
potential should be done independently. Generalization of our approach to this case should 
be possible and we plan to do it in the future. 

Another important feature of our theory is that it is gives, by construction, finite cor- 
relators for the vertex operators, corresponding to the fermions in the quantum mechanical 
picture. This property, natural for the effective theory of non-relativistic fermions, looks 
strange from field-theoretical point of view. It cannot be present in the CFT. In fact, our 
theory is not conformal because of the presence of an unusual background term. Modified 
Liouville field in our theory is in some sense auxiliary. If one integrates it out, one gets an 
effective action for the field X. This action is complicated because it should give UV finite 
correlators of the vertex operators. If one uses the simplest approximation for the integral 
over the Liouville field (as we discussed at the end of Section m| . one gets just a sharp wall - 
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Dirichlet boundary state making X chiral. However, this approximation is non-homogeneous 
in z-plane, as it was discussed above (Section |2I). Indeed, near the boundary of the droplet 
(Fermi sea) the step function, which is fermion density, corresponding to the naive large N 
approximation, is always far from the real density profile (in a sense that for any N there 
is z such that the difference between the exact pn and the step function is of order unity). 
If one uses more careful approximation for the Liouville field, the effective action of the X 
becomes more complicated after we integrate out modified Liouville field. 

Thus what we have obtained a UV-finite theory, describing bosonization of free fermions. 
Because of its apparent similarity with the two-dimensional string theory, it is natural to 
ask: "is this some sort of string theory?" . The answer is the following: if there is a string 
theory description of our theory ()4.17|) . it is a string theory in its target space formulation, 
a string field theory. Then it is interesting to understand what should be the world-sheet 
description of such a theory, if it exists. It would be also interesting to see how such a 
description will look like for the inverse harmonic potential. As we discussed above, we 
expect the derivation to be very similar and the answer to differ from the case at hand by 
real instead of imaginary Liouville field and Lorentzian signature. Other features of our 
theory should remain unchanged. It should be possible to compare the action which appears 
after integrating out the Liouville field with the Das-Jevicki action. 

Acknowledgments 

We acknowledge useful communications with A. Abanov, S. Alexandrov, V. Kazakov, I. Ros- 
tov, E. Martinec, N. Nekrasov, J. McGreevy, H. Verlinde, S. Wadia, P. Wiegmann, A. Zabrodin 

A Careful treatment of the ultra-violet and infra-red 
regulators 

In this section we carefully treat both IR and UV regularizations of the theory (j2.12j) . For 
the grand canonical ensemble of the Coulomb plasma (j2.9|) 



we rewrite the Hamiltonian (j2.4j) in terms of a correlator in a free boson field theory: 



^"^When this paper was finished, the work |33| appeared, which suggested a string theory, which can serve 
as a possible candidate for the world-sheet description of our theory. A remarkable feature of the theory 
is that its world-sheet action is quite similar to that of H4.17II . The relation between the two theories is 
currently under investigation. 




(A.l) 




(A.2) 
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where now the average is taken with respect to the Euchdean Gaussian action of a free 
massive boson (thus introducing infra-red regulator into the theory): 



Sri 



27ip 



(A.3) 



with a short-distance cutoff A (see Eq. ()A.5|) below) implied. Formula ()A.2|) is exact and 
follows from the the Gaussian integration formula 



N 

n 

\j=i 



exp 



1 ^ 



(A.4) 



where the correlation function of the field d> is 



{<Piz,z)m)s. 



2/3Ko{m\z\), 



2/3 In f A 



A\z\ > 1 
A\z\ < 1 



(A.5) 



Here -^'0(2;) is the MacDonald function with asymptotics given in Eq. (HH). 

- In X — 7 + In 2 + 0(x), x —>■ 0, 



Kn(x) 



^e--[l+o(l/x)], 



X ^ 00 



(A.6) 



For the regularized version of the partition function Ziy{P), eq. ()2.3p . the representation 
2)) gives 



N 



2^ / TTgj0{2j,2j)-/3|^j|2 



(A.7) 



From the asymptotic ()A.6|) one finds that 

Tat = lim 



m->0 



^N,reg 



(A.8) 



where tn is given by ()2.2j) . A reason for having the infinitesimal mass m in the action ()A.3 
is clear from eq. ()A.4jl . For exactly zero mass the integration over the zero mode of field ( 
would make the right hand side exactly zero. 

Substituting the representation ()A.7|) in ()A.1|) we find 



Z 



Zs ^ N\\m 



f3N 



(A.9) 



The sum over N is easily performed yielding 



1 



Za 



(A.IO) 



23 



where 



iL,re,M = 7^ fd'z \d<f)d<j) + m^^ - /xe^^-^l^ll (A.ll) 
iTrlJ J L J 

Equations ()A.10|) and ()A.11|) define the path integral representation for Z. The couphng 
constant 

A\^ 



/x= - 27r/?e (A. 12) 

\m J 

has the dimension of inverse length squared. It is a running coupling constant satisfying the 
RG equation 

(ilogA ^ ' 

We are interested in the m — limit of the grand partition function ()A.1|) . In this limit 
grand partition function ()A.1|) possesses a remarkable property: grand canonical ensemble 
gives an exact approximation for the averages calculated in the microcanonical ensemble. 

Indeed, assuming iV ^ 1 (where iV is a fixed average number of particles) one can 
perform the sum in eq. ()A.9|) in the saddle point approximation. The saddle point equation 
defines the chemical potential x = log'C terms of the average particle density N 

X = ^^N-N\n + log iV (A. 14) 

where f^Ar is a thermodynamic potential of the Coulomb plasma ()2.H|1 . For iV ^ 1 it has the 
following asymptotic (see e.g. [^): 

= -]^N'^\iyN ^'^^N'^ ^W[N) (A.15) 

This can be understood by computing the free energy of the uniformly charged disk of radius 
in the uniform background of the unit charge (the subleading term W{N^ = 0{N \n N) 
is the correlation energy). The fluctuation of the number of particles 6N around the saddle 
point is given by the standard thermodynamics formula 

iSNf ^ if ,A.16) 

Substituting in this equation the result eq. ()A.14|) for the chemical potential x and using 
formula ()A.15|) we find 

This is a remarkable formula. It suggests that in the limit m — > the number of particles in 
the plasma does not fluctuate and the grand canonical ensemble gives an exact approximation 
for the averages calculated in the microcanonical ensemble. 

Last in this section we compare the average calculated in canonical and the grand canoni- 
cal ensemble. Consider an observable 0{zi, . . . , zn)- Denote its average in the microcanonical 
ensemble 

= 7 ^ / d^^i . . . d'zN e-''^--0(zi, ...,zn) (A.18) 

^N,reg[P) J 
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The average in the Grand canonical ensemble will then be given by 

N=0 

Performing the sum in the saddle point approximation we find 

(O) =Or, + mf\-^,0^. (A.20) 
By virtue of Eq. ()A.17j) the correction vanishes in the limit m : 

(C) ^ Ojv, m ^ (A.21) 

B Expansion of the density 

In this section it will introduce an additional parameter h = 1/N to keep the area of the 
droplet of the ground state finite and keep (3 = 1 for simplicity. The asymptotic expression 
for the one-particle density matrix is obtained from the exact formula ()H.2H) 

Using the small expansion for the logarithm of the r-function 

log K7^^N{tk,tk) = log Cn + ^y^^^kh + o{\tkf). (B.2) 



fc>i 



(where Cat = tviP'tt)' r-function in ()B.1|) we find 

1 / . . \ C 



It is instructive to compare the asymptotic formula ()B.3|) with the exact result in the P = 1 
case, which can be found from comparing the expression ()B.1|) with the exact density matrix 
of free fermions 



p(U) = e-i^^^'Y.^, (B.4) 

A;=0 



The resulting expression for the r-function reads 

1^ fcc*' K'')-^N2^(fq_jA\^kik y^-^i 



One can see that for h = 1/N the coefficients of the sum ()B.5jl are rapidly decaying with k. 
Simple asymptotic analysis shows that for small k the Stirling formula can be used for all 
factorials giving after some algebra 



exp 



{N-k)\ 



e k^ 

' 1 ■ 2Ar2 + 2 ■ 3iV3 ^ • • • 
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(B.6) 



One can see that this expression is exponentially small for p > \/N and therefore only the 
leading term in the exponential can be kept. This gives 

^r„(-^.-^)=.C.ge-i^^ (B,7) 

Comparing equation ()B.7|) with ()B.3|) one can see that the asymptotic expansion ()B.3|1 breaks 
down for p ~ \fN . The same is true for the small expansion ()B.2jl . 
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